Defect formation in the Swift-Hohenberg equation 
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We study numerically and analytically the dynamics of defect formation during a finite-time 
quench of the two dimensional Swift-Hohenberg (SH) model of Rayleigh-Benard convection. We 
find that the Kibble-Zurek picture of defect formation can be applied to describe the density of 
defects produced during the quench. Our study reveals the relevance of two factors: the effect of 
local variations of the striped patterns within defect-free domains and the presence of both point-like 
and extended defects. Taking into account these two aspects we are able to identify the characteristic 
length scale selected during the quench and to relate it to the density of defects. We discuss possible 
consequences of our study for the analysis of the coarsening process of the SH model. 

PACS numbers: 64.60.Cn, 47.20.Bp, 47.54.+r, 05.45.-a 



The formation of topological defects in symmetry- 
breaking phase transitions is a very generic phenomenon 
in physics, and can be studied analytically and exper- 
imentally in different condensed matter systems [jj, 
An example is the onset and formation of stripe pat- 
terns in Rayleigh-Benard convection 0, • In this pa- 
per, we focus on the Swift-Hohenberg (SH) model '3| 
for this process. Once above the convective threshold, 
the system develops a labyrinthine morphology, consist- 
ing of domains of stripes which are oriented along arbi- 
trary directions Between those domains, the system 
displays several types of topological defects such as grain 
boundaries, disclinations and dislocations. This struc- 
ture orders with time basically by grain boundary relax- 
ation and defect annihilation. Similar to other models, 
in which this coarsening process is self-similar 0, any 
linear scale of the structure is expected to grow as a 
power law in time ^ ~ t^^^. However, simulations of 
sudden quenches of the SH equation (7ii. .8. .9. ,10 ^ iT H Il^ 
have revealed that the observed exponent z is sensitive 
to non-universal model features such as the quench depth 
or noise strength. Moreover, different definitions of the 
length scale have led to different exponents with values 
reported in the interval 2 < z < 5. This apparent ab- 
sence of self-similarity in the coarsening is also found in 
related experiments of electroconvection and diblock 
copol yme rs [T^ . Multiscaling and/or defect pin- 

ning [Tl| have been proposed as being responsible for 
the scattered value of z, but no general picture has been 
reached so far about the true nature of the coarsening 
process. 

In this paper we consider a complementary, but re- 
lated aspect of the non-equilibrium dynamics of the SH 
equation. Namely, we are interested in the formation 
of defects in a finite-time quench (annealing). Interest- 
ingly, some features of finite-time quenches were consid- 



ered in some early works comparing the SH model with 
Rayleigh-Benard experiments . Specifically, we study 
situations in which the control parameter, the reduced 
Rayleigh number e = {R — Rc)/ Rc, is swept smoothly 
over the bifurcation point, e = 0. Our study is thus 
related to earlier papers on defect formation in nonequi- 
librium second-order phase transitions 0, HtL IT^ . The 
theoretical picture which is believed to be applicable in 
this case is the Kibble-Zurek mechanism because of 
critical slowing down close to the transition point, the 
system can not follow the external change of the critical 
parameter and the dynamics are suspended in an inter- 
val [— e, i] around the transition. After this suspension 
a characteristic length ^ is selected and is found to scale 
with the rate of increase, ^, of the critical parameter 
like ^ ~ A^"'''- In particular, mean-field theory predicts 
e ~ /i^/^ and ^ /i^^/^ for models of the Ginzburg- 
Landau type with 0{N) symmetry 16]. This length ^ 
sets an initial density of defects, p, directly after the 
quench. Within a Gaussian approximation one expects 
p ^ for point-like defects and p ^ for line defects 
in two dimensions Thus, the Kibble-Zurek mecha- 
nism predicts that the density of defects observed directly 
after the quench scales with p as well, as confirmed in 
various models [l7lllll20l|. 

The purpose of this paper is to confirm the validity 
of the Kibble-Zurek scenario for the SH model. To this 
end, we identify the characteristic length ^ selected dur- 
ing the annealing process and study its relation to the 
density of defects. We model the annealing protocol of 
the Rayleigh-Benard system by the SH equation in two 
dimensions and in dimensionless variables 0, Q 
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where e{t) — pt. This corresponds to an experimen- 
tal situation where the temperature difference between 
the upper and the lower plate of the convection cell is 
increased linearly in time. The order parameter field 
(/)(x, t) is related to the vertical fluid velocity. The 
last term in Eq. Q is a stochastic forcing term, with 
{rj{x,t)r]{x',t')) = 2F6{x - x')S{t - t'), where the noise 
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FIG. 1: Variance of the order parameter, ((^■^), obtained 
from the simulations (sohd hnes) compared with the hnear ap- 
proximation given by Eq. Q (dashed hnes) with = 10~^, 
^ — 10~^ and /i — 10~^ from left to right. Inset: numer- 
ical values of e as a function of fi obtained implicitly from 
{(j)^) — Si with 5 = 1/2. Solid line is the prediction of Eq. ||SJ. 



strength F = 5 x 10~^^ is compatible with typical exper- 
imental values 0,^3 • Our simulations are performed on 
a square lattice of 512 x 512 nodes with lattice spacing 
Ax = 7r/4, corresponding to 8 lattice points per ideal 
wavelength {qo = 1). Initial conditions are (/)(x, to) = 
at the instant tg given by e{to) — —1/2. 

As the critical parameter e is increased in time the 
magnitude of the order parameter remains close to = 
until well after the onset of the instability at e = (see 
Fig.^. At some later instant £ = £ > the field abruptly 
jumps towards its symmetry-broken quasi-equilibrium at 
(0^) = 2e{t)/3 and spatially periodic modulations of the 
signal are created, see Fig.|21 We identify the instant e at 
which dynamics are resumed as the time when (0^) = 6e 
p8] [with ,5 = 0(1), meaning that the first and third 
term on the right-hand side of Eq. Q are of equal im- 
portance]. Results are presented in Fig. ^ where we find 
£ ~ ^o.48±o.oi^ which compares well with the Kibble- 
Zurek predicted scaling e ~ fi^/^. It is interesting to note 
that the same scaling was found in ramp experiments of 
Rayleigh-Benard convection 15j , where the temperature 
difference between the upper and lower plate of the con- 
vection cell was increased linearly in time. This setup is 
similar to the annealing protocol in our simulations sug- 
gesting that the scaling found in those experiments can 
be explained within the Kibble-Zurek scenario as well. 

Since the order parameter remains small up to £ = 
£, the observed scaling can be obtained from a simple 
linear approximation of Eq. Fourier transforming 

the linearized SH-equation yields 

dt^t) = [^,t - {ql - q2)2] 0q(t) + f,^{t). (2) 

where 4>q(t) is the Fourier mode of the order parame- 
ter field with wave vector q. Thus, the structure factor 
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FIG. 2: Typical configurations of the order parameter field 
for different values of e. We depict an area of 256 x 256 lattice 
points. Upper row: jU = 10~^, lower row fi — 10~^. Black 
(white) points correspond to (p > {(j) < 0) . Rightmost pan- 
els display the defect structures obtained from the amplitude 
signal (i? — 0.25, see text). 



Siq,t) = (|</)q(t)P) is given by 



S'(g,t) = 2Fe2'^'(*) f 

J to 



(3) 



where ujq{t) = ^^j.{t'^ -t^) - {q^ - q^)'^{t-tQ). The struc- 
ture factor S{q, t) is peaked around q ~ go (see Fig. |3|d) 
and the time dependent width of this peak, T(t), is ob- 
tained from S[qo ± ^T{t),t] = ^S{qo,t). We find that 
r(t) satisfies the implicit equation 
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where a = 1 -I- 0{T^). Assuming that S{q,t) is sharply 
peaked around qo (i.e. F ^ 1) and that it can be approx- 
imated by a squared Lorentzian around go 0.0 (see Fig. 
ISJs) we write 



qaS{qQ,t) 
8(V2- 1)1/2 



(5) 



Combining this result with the implicit equation 
6s, we find that e satisfies 



fiH{2F)-'6qoCo. 



(6) 



where C'o ~ 4.93 is a numerical constant. Thus, the linear 
approximation leads, up to logarithmic corrections, to 
the scaling behavior e ~ /i^/^ predicted by the Kibble- 
Zurek scenario and confirmed by our simulations (see Fig. 

m 

We now proceed to identify the typical length scale 
selected by the dynamics during the quench. A canonical 
measure for this leng th is the width of the structure factor 
11 H H E in Ell : in a defect-free domain the order 
parameter field takes the form 

0(x, t) ~ ^(x, t) cos(q(x, t) ■ x)-|-higher harmonics (7) 
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FIG. 3: a) Width of the structure factor, F, after the quench 
compared with the local variations of the wave-number Aq. 
Solid line is the approximation of Eq. (|HJ , dashed lines are fits 
to power laws F(£ = i) ~ /■2^±o «^and F(e = 1) ~ / iBio w. 
b) Collapsing plots of S{q) as a function of the rescaled vari- 
able X — ((f — go)/(QoF) at e = e for /i = 10^^ (dotted line), 
= 10""^ (dashed line) and = 10""* (long dashed line); the 
solid line is a squared Lorentzian curve 1/(1 -|- x^Y . 

where A(x, t) and q(x, t) vary slowly within the domain. 
Thus, in the absence of any defects one expects the struc- 
ture factor S{q, t) = (l'^q(i)P) of an infinite system to be 
a delta function centered around q = Qq- When defects 
are present, the structure factor broadens around qo and 
the width F at half-height can be used as a proxy for the 
length, i.e. we define ^ = F^^. Fig.|3Ji shows the data for 
F at e = e and s = 1 numerically obtained from a circu- 
lar average of the structure factor. We find that S{q, t) 
has a scaling form S{q,t)/ S{qo,t) ~ /[(q^ _ q^)/{qQT)] 
(see Fig.Eb) at least for a given (scaling) interval around 
g = go — 1- Results are thus independent of the thresh- 
old (half-height) used to obtain F, provided the analysis 
is performed within this scaling region. At e = s we find 
^ ^-o.24±o.oi again confirms the Kibble-Zurek 

prediction ^ ^ Note that this scaling can be re- 

covered within the linear approximation by combining 
Eqs. Q and (O to obtain 



m ^ M^^y || [H5qoCos'/'{2F)-')]-y^ (8) 

which agrees with the Kibble-Zurek result up to logarith- 
mic corrections and compares quantitatively well with 
the numerical data (see Fig. Eli). 

We find, however, that right after the jump the ab- 
solute value of ^ decreases and then rapidly saturates 
to a value smaller than the one at the jump. The 
data at any later instant (e.g. e = 1) deviates from the 
Kibble-Zurek prediction in the regime of small /U, in fact 
in the interval 10"'' < /i < 10"'^ it can be fitted to 
^ ^ F^^ ^ ^-o.iSito.oi^ understand this apparent 
inconsistency, we first note that coarsening of the de- 
fect structure cannot be responsible for this effect, since 



£,{e = e) > £(g = 1). Instead, a possible explanation is 
given in 0, 0|: the broadening of the spectrum and 
hence the observed value of F is not only related to the 
spatial distribution of defects, but also to variations of 
the local wavenumber g(x, t) within the defect-free do- 
mains. Using the numerical procedure suggested in [2ll | 
we find that although the distribution of local wave- 
numbers q(x, t) in Q is centered around the expected 
value ^(x, t) = qo, it exhibits a finite width indicating 
variations in the actual values of q(x, i). To quantify 
this effect we measure the RMS of those fiuctuations, 
Aq = ([g(x, t) — qo]"^)^^^, where the average extends over 
defect-free domains only. We find that Aq is a significant 
fraction of F for small values of /i, where F is small (see 
Fig. Ot)- Although in principle the observed value of F 
can be a complicated convolution of both the local varia- 
tions of the value of q(x, t) and the presence of defects, we 
have tried a simple linear Ansatz to disentangle both con- 
tributions, i.e. F = Fdcf + Aq, where we assume that Fdcf 
is directly related to the spatial distribution of defects. In 
fact, the data for Fdef at e = 1 obtained using this simple 
Ansatz can be fitted to a power law F^of ^ ^-o.24±o.oi 
and rescaled by a constant factor to collapse well with F 
taken at £ = £ as shown in Fig. This indicates that 
the Kibble-Zurek scaling f ~ might still hold at 

£ = 1, although it appears to be masked by local varia- 
tions of the wavenumber 22]. 

An independent way of checking the Kibble-Zurek pre- 
dictions is to measure the density of defects, p, created 
during the quench. Several methods have been devised 
to identify defects numerically _8, 10^ 0, Here 
we present results for two independent measures of the 
number of defects: the first one is based on the local 
amplitude of the order parameter field , while the sec- 
ond one relies on the local curvature of the observed 
stripes [Tll |. Far away from any defect the order pa- 
rameter field 0(x, t) is of the sinusoidal form ||7J), so 
that the local amplitude can be numerically estimated 
as A'^(x,t) ~ 0(x, t)^ -I- (V(?!)(x, t))^/gQ. In equilibrium, 
it is found that A^(x, i) ~ Aq = Ae/Z to lowest or- 
der in £ for a set of parallel stripes The density 
of defects, pa, can thus be calculated using the filter 
|[A^(x, i) - Aq]/Ao| > -d, where d is some threshold. 
The second method consists in numerically extracting 
the probability distribution P(k, t) of the local curvature 
k(x, t) — |V • n(x, t)| of the stripe patterns 0. Here 
n — V(/)/|V0| is the unit vector normal to the lines of 
constant (p. For a defect-free pattern, the curvature is 
zero, so again defect points can be identified as those 
with k(x, i) > 0, where 8 is a given threshold. Thus 
the density of defects is pc — Jq P{K,t)dK. In Fig. El we 
display typical configurations of the SH signal and the 
corresponding defect structures obtained using the first 
filter, demonstrating that this method can detect the de- 
fect structure efficiently. The thresholds and are 
carefully chosen to be within the finite interval in which 
our results for pa and pc are (up to a constant prefactor) 
independent of the particular choice of and Q. 
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FIG. 4: (a) Density of defects pa and pc (rescaled to collapse 
with Pa) obtained at e = 1 (the data for fi > 10~^ is taken at 
e — 2), and Fdof = r(e = i) (in arbitrary units) as a function 
of fj,,. Solid line corresponds to p ~ /jP'*^ and dashed line 
is Fdot ~ /i"'^*. (b) Density of defects as a function of Fdof- 
Dashed line is p ~ ^dei- 

As shown in Fig. 0] the numerical values for pc can be 
rescaled to agree very well with the density of defects 
obtained from the amplitude signal, pa, so that (up to 
a constant factor) both methods independently yield the 
same density of defects. For relatively large values of p,, 
we find that the density scales like p ~ /x^^^ and also that 
p ^ r^^j, as shown in Fig.^ Consequently, if Fdof is re- 
lated to the density of defects, those should be point-like 
for quick quenches. A simple inspection of the configu- 
rations after the jump confirms this picture, see Fig. |21 
we observe that for small values of p, most of the defects 
are point-like, while only a minor fraction is spatially ex- 
tended. Only at small values of the sweep rate does the 
density of defects deviate from the p ~ ^'^^t /^^^^ scal- 
ing, and we observe an increasingly significant amount 



of extended defects. If extended defects dominate over 
point-like defects, a crossover to p ^ Fdcf ~ M^^'' should 
be expected for small p. Although a departure from the 
p ~ p^/^ behaviour is evident for slow quenches, our com- 
putational facilities at present do not allow us to give an 
accurate confirmation of this possible crossover. 

In conclusion, wc have confirmed the validity of the 
Kibble-Zurek picture for the creation of defects in the 
annealing of the SH equation. Specifically, we have iden- 
tified a length scale ^ ~ FJ^^ directly related to the den- 
sity of defects, and found that it scales with the sweep 
rate as ^ /i"^/"*. For large values of p, where most of 
the created defects are point-like, we expect p ~ a 
result which is observed in our simulations. However for 
slow quenches, where an increasing fraction of extended 
defects is produced, we expect a deviation from this scal- 
ing law and a possible crossover to p ~ . We hope this 
work will stimulate other simulations of the SH model in 
order to confirm this picture. 

Finally, we note that in studies of the coarsening in the 
SH model mil , it is observed that at large times defects 
are extended and that thus p ^ £,~^- Our study suggests 
that the presence of point-like defects might be relevant, 
and thus some care has to be taken in the SH model when 
comparing the density of defects with the characteristic 
length observed in the system. We hope that taking into 
account this aspect together with the identification of 
the length scale associated with the density of defects 
as suggested in this paper might be useful to shed some 
light on the question of self-similarity in the coarsening 
process after an sudden quench of the SH system . 
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